Inner and outer edge states in graphene rings: A numerical investigation 
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We numerically investigate quantum rings in graphene and find that their electronic properties 
may be strongly influenced by the geometry, the edge symmetries and the structure of the corners. 
Energy spectra are calculated for different geometries (triangular, hexagonal and rhombus-shaped 
graphene rings) and edge terminations (zigzag, armchair, as well as the disordered edge of a round 
geometry). The states localized at the inner edges of the graphene rings describe different evolution 
as a function of magnetic field when compared to those localized at the outer edges. We show 
that these different evolutions are the reason for the formation of sub-bands of edge states energy 
levels, separated by gaps (anticrossings). It is evident from mapping the charge densities that the 
anticrossings occur due to the coupling between inner and outer edge states. 

PACS numbers: 73.43.-f, 73.23.-b, 73.63.-b 



I. INTRODUCTION 

Graphene is a bona fide two-dimensional material 
showing a great versatility due to its unconventional elec- 
tronic properties and promising applications to nanoelec- 
tronics [TJ |2]. Among the promises is the possibility of 
structuring graphene at a mesoscopic length. Indeed, 
some groups have already demonstrated that graphene 
can be cut in many diS^erent shapes and sizes, opening the 
door to the fabrication of graphene nano-devices through 
the impressive experimental obtention of graphene quan- 
tum dots [31 m [5], quantum rings [6] and even antidot 
arrays [7]. This perspective leads to interesting scenar- 
ios since the electronic properties of graphene are deeply 
influenced by its size and shape: as it is well known for 
over 10 years, graphene nanoribbons have different prop- 
erties depending on their width or on their edge termi- 
nations [8j m Uni HI]. Besides the nanoribbons, some 
theoretical works have also addressed the effects of con- 
finement on the electronic structure of graphene quantum 
dots (flakes) with different geometries, sizes and types of 
edge [12j [131 E] • Recently, the energy levels of graphene 
quantum rings [15l [16] and of graphene antidot lattices 
[Tfl [TB] have also been theoretically investigated. 

In this paper we numerically analyze the energy spec- 
tra as a function of magnetic fleld {B) of graphene quan- 
tum dots and graphene quantum rings, with focus on the 
complex evolution of edge states in the graphene rings. 
Here we explore the effects of the interplay among dif- 
ferent degrees of freedom given by size, geometry and 
edge symmetries on the electronic properties of these 
graphene nanostructures. We consider quantum dots and 
rings with six (hexagons), three (triangles) and two- fold 
(rhombus-shaped) rotational symmetry, with zigzag or 
armchair edges. Afterwards we consider round dots and 
rings, whose edges are not so simply: they are cut in a 
way to approach the circular geometry. Our attention is 
concentrated on the continuum limit [H] of the energy 
spectra as a function of the magnetic flux of these struc- 



tures. Edge states appear with energies between consec- 
utive Landau levels (LLs) in such spectra, as could be 
initially expected [20 ,. However, the interplay between 
two different edges showing distinct local structures (the 
quantum rings can be seen as graphene structures con- 
taining an antidot, which introduces an inner edge to the 
system) leads to some surprising subtleties. We observe 
that the presence of the antidot introduces additional 
edge states, with a different evolution under B: their en- 
ergies are increased with increasing B. For a better un- 
derstanding of this behavior, the electronic densities of 
these states are mapped, and we show that edge states 
that rise in energy with B are located in the internal 
edges of the ring structure. In this way, we show that in- 
ner and outer edge states give origin to the formation of 
sub-bands separated by energy gaps in the region of the 
spectra between LLs. The anticrossing of levels, which 
defines the sub-bands, occurs due to the inter-edge cou- 
pling of states. The formation of sub-bands is highly 
influenced by symmetry properties, and also by size ef- 
fects, i.e., the relation between the ring width and the 
magnetic length. 

As will be seen throughout the paper, the choice of 
quantum rings is strategic since the states within the edge 
states sub-bands can be perfectly associated to either in- 
ner or outer edges, or to a coupling of both edges of the 
ring structure, therefore enabling a good framework for 
studying the influence of the edges and edge junctions on 
the electronic structure and charge distribution. 



II. MODEL 

We use a tight-binding model for a finite two- 
dimensional honeycomb lattice, considering nearest 
neighbors hoppings [HI [21]. The following non- 
interacting Hamiltonian is considered: 
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H = Y^ e4c, + t (e'^-cfc, + e-^^-cJc,) (1) 

i <ij> 

where Ci is the fermionic operator on site i. The per- 
pendicular apphed magnetic field is included by means of 
Peierls substitution, which means a complex phase in the 
hopping parameter {t=2.7 eV): = 2iT{e/h) A-dl . 
In the Landau gauge, where the electromagnetic vec- 
tor potential is defined as A = {0,Bx,0), one obtains 
— along the x direction and (f)ij — ±7r(a:/a)<i>/<i>o 
along the direction. The magnetic flux (<I>) per mag- 
netic flux quantum (i>o — h/e) is defined as: ^/^o — 
Ba'^V3e/{2h), and we use a=2.46A as the lattice con- 
stant for graphene. The on-site energies are taken as 
e, = 0. 

To consider the ring geometries, a central region of ab- 
sent atoms (antidot) is defined in the structure by setting 
the hopping parameters to zero for the absent atoms and 
the on-site energies at the position of these atoms equal 
to a large value outside the energy range of the spectra. 
The magnetic field we consider is not limited to the cen- 
tral region of the ring, but is homogeneously applied to 
the entire structure. By exact numerical diagonalization 
of the Hamiltonian, the energy spectrum as a function 
of magnetic fiux is calculated for different geometries of 
quantum rings. 

III. EDGE STATES OF GRAPHENE RINGS 
A. Antidot effects on the energy spectrum 

We start by calculating the energy spectrum of the 
structure shown in Fig. 1(a): a finite hexagonal lattice 
forming an equilateral triangle with zigzag edges. Call- 
ing Nout the number of individual hexagonal plaquetes 
along each side of the triangle, the total number of car- 
bon atoms in this structure is N^^t + ^^out + 1 [H] • The 
energy spectrum for such a triangular graphene quantum 
dot with iVo„t=45 is plotted in Fig. 1(c) as a function 
of magnetic flux. One can clearly observe the forma- 
tion of the low energy LLs: the n=0 LL at zero-energy, 
and the n=+l and n=+2, with their square-root depen- 
dence on magnetic field, typical from graphene systems 
(the spectrum is symmetrical with respect to the zero- 
energy). Also, one can see the expected presence of edge 
states between consecutive LLs and observe the evolution 
of these edge states with magnetic fiux until coalescing to 
the LLs [ID] . The side length of the triangular structure 
is simply given by aNout, where a=2.46A is the lattice 
constant. So, for the case considered here of iVo«t=45, 
the side length of the triangular dot is «llnm. 

We then take this nanostructure as a starting point to 
develop a triangular quantum ring just piercing a trian- 
gular hole (antidot) in the middle of it. The inner edges 
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FIG. 1: (a) Triangular graphene quantum dot with zigzag 
edges, with Nout=^^- (b) Triangular graphene quantum ring, 
with 7Vout=45 and Nin=12. (c) Energy spectrum as function 
of the magnetic flux for the structure in (a), (d) Energy 
spectrum as function of the magnetic flux for the structure in 
(b) 



of this ring are also zigzag. To define the size of the an- 
tidot, we call the number of hexagons at each side 
of the internal removed triangle. The total number of 
atoms in this quantum ring is now Nl^^ + A{Nout — Nin) — 
Nf^ + QNin- In Fig. 1(b) there is a representation of such 
a ring for Nout=^^ and 7Vi„=12, and the corresponding 
energy spectrum is shown in Fig. 1(d). The interesting 
observation is that the presence of the antidot gives ori- 
gin to additional edge states with a different evolution 
with magnetic field: states that go up in energy as the 
magnetic fiux is increased. Indeed, looking at vacancies 
in graphene |T9], the locahzed states around such defects 
also rises in energy with increasing magnetic field, since 
vacancies are actually minimal antidots. It is also clear 
that with the introduction of the antidot, the formation 



of the n >0 LLs starts at higher magnetic fields when 
compared to Fig. 1(c), due to the involved interplay be- 
tween the inner and outer edge states. On the other 
hand, this interplay seems to anticipate the formation of 
the central LL, in these peculiar zero magnetic field limit 
showing edge states at the Dirac point due to the zigzag 
structure of the edge [S] . 

To analyze in more details how the energy levels of the 
edge states evolve with magnetic field, in Fig. 2(a) we 
zoom in the energy scale of Fig. 1(d). It now becomes 
evident the formation of edge states sub-bands, separated 
by energy gaps that get smaller with increasing field. One 
can also note that each of these sub-bands contains three 
crossing energy levels for this triangular graphcnc ring. 



B. Different evolutions for inner and outer edge 
states 

In order to gain a deeper understanding of this quite 
complex evolution of edge states, in Figs. 2(b-d) we look 
to the wave functions amplitudes of specific edge states. 
The arrow (b) in the spectrum in Fig. 2(a) is pointing 
to an edge state whose energy is reduced with increasing 
B. This state is mapped in Fig. 2(b), and clearly is 
localized at the outer edge of the triangular ring. The 
radii of the circles plotted are directly proportional to the 
wave function amplitude on each site, and we can observe 
a symmetrical and quasi homogeneous distribution over 
the edges, with higher concentrations at the outer most 
lattice sites (and always on the same sublattice, in this 
case of zigzag edges). A high charge density accumulation 
is also observed close to each corner forming part of a 
second charge density belt. 

The arrow (c) in Fig. 2(a) points to one of those states 
that go up in energy with B, whose wave function is 
mapped in Fig. 2(c). In agreement with the observation 
that those states going up in energy appear only when 
the antidot in considered in the lattice, this is an edge 
state clearly located at the inner edge of the triangular 
quantum ring. Its electronic charge density is homoge- 
neously distributed over the innermost lattice sites with 
decreasing amplitudes when approaching the corners. 

This system has a three-fold rotational symmetry. The 
outermost atoms at the outer edges are all from the same 
sublattice except for the three atoms located at each cor- 
ner. The atoms at the innermost edges are all from the 
same sublattice, including that located at the corners, 
and this is an important difference between the inner and 
outer edges in the triangular zigzag quantum ring. The 
charge density in the innermost edge is also sublattice 
modulated but with the charge density dominantly on a 
different sublattice than at the outermost edge. 



3 



0.06 



0.04 

0.02 




FIG. 2: (color online) (a) Zoom in the energy scale of the 
spectrum shown in Fig. 1(d), showing now only the first few 
low-energy states and their evolution with magnetic flux, (b- 
d) Electronic charge distribution of the selected edge states 
indicated by the arrows and corresponding letters in the spec- 
trum. Three typical behaviors are clearly defined: (b) an 
state whose energy is reduced with B is an outer edge state, 
(c) an state whose energy is increased with B is an inner 
edge state, (d) at the anticrossing levels the wave functions 
are distributed between the inner and the outer edges, indi- 
cating a coupling between both edges. The radii of the circles 
are proportional to the amphtude of the charge density. 

C. Coupling between inner and outer edge states 

An interesting observation emerges from mapping the 
charge density of a state situated at an anticrossing, like 
the state indicated in the spectrum by the arrow (d) and 
with the charge density depicted in Fig. 2(d). It can 
be observed that the wave function has amplitudes con- 
centrated on both the inner and the outer edges of the 
ring. This indicates that a coupling between edge states 
from the inner edge and from the outer edge is taking 
place. Consistently to this picture of inner and outer 
edge states seeing each other and getting coupled, we 
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observe that the higher the magnetic flux, the smaller 
are the energy gaps between the sub-bands. This can be 
attributed to the fact that increasing the magnetic flux 
<i>/<i>0) the magnetic length Ib of the states gets smaller, 
reducing the chances of coupling. 

For a more quantitative comparison, the width (dis- 
tance between outer and inner edges) of the triangu- 
lar ring we are considering (with iVo„t=45 and iVi„=12) 
is 24.1 A and the magnetic length is determined by: 
Ib = ^Jh/eB = 0.913A/^$/$o- In this way, for a 
flux <i>/<i>o=0.02, for which there are no energy gaps 
in the scale observed in the spectrum of Fig. 2(a), 
we have Ib=6.46A. Reducing the flux, for example for 
<I>/$o=0-01, where energy gaps already start to appear, 
the magnetic length is Zb=9.13A. For $/$o=0.005, a re- 
gion of flux where the energy gaps are more clearly de- 
fined, we have Zb=12.9A, a value corresponding to ap- 
proximately half the width of the triangular ring, and so 
compatible with the suggested coupling between outer 
and inner edge states. 

We recall the fact that we are showing typical charge 
density plots: any other chosen edge state that goes down 
(or up) in energy has a very similar charge density dis- 
tribution to that shown in Fig. 2(b) (or 2(c)), while 
any state at an anticrossing shows wave function concen- 
tration on both inner and outer edges, similarly to the 
distribution observed in Fig. 2(d). It is interesting to no- 
tice that the coupling of both edges does not break the 
sublattice modulation of the charge density at each edge, 
but the state as a whole is now sublattice mixed. 



IV. WIDTHS, EDGES AND CORNER EFFECTS 
ON THE QUANTUM RINGS 

A. Widths and sub-band gaps: tuning the inner and 
outer edge coupling 

We now turn our attention to hexagonal graphene 
quantum rings, first to compare the energy spectrum of 
this other geometry with the one from the triangular ring 
shown previously, and second to show the interesting ef- 
fects of varying the width of the quantum ring on the 
formation of the edge states energy sub-bands. In Figs. 
3(a) and 3(b) there are representations of two hexagonal 
quantum rings, with different widths. We once again con- 
sider a geometry with all zigzag edge terminations. The 
total number of atoms in an hexagonal zigzag graphene 
ring like these is 6iV^„j — 67Vf„, where Nout and 7Vi„ are 
the number of hexagonal plaquetes along each side of the 
hexagon and removed hexagon, respectively. Again, the 
total length of each side of the structure is just given by 
the number Nout or Nin times the lattice constant a. For 
the ring in Fig. 3(a) we consider A^o„t=21 and Nin—7, 
while the ring in Fig. 3(b) has Afo„t=21 and Nin=12. 

The energy spectra as a function of magnetic flux of 



these two rings are shown in Figs. 3(c) and 3(d), respec- 
tively. Comparing these spectra with one of an hexagonal 
quantum dot (without the antidot in the middle) [13], it 
is evident that the ring geometry introduces energy sub- 
bands separated by energy gaps, exactly as in the case 
of the triangular quantum ring. However, it can now be 
observed that each band of the hexagonal structures has 
six energy levels instead of the three levels observed in 
the triangular structures. We note that this follows the 
rotational symmetry fold number of the hexagonal ring 
structure. Here the sublattice of the outer and inner most 
edges alternates from one sublattice to the other going 
from one arm of the hexagon to the next |13} 115] . 

The width of the ring (distance between outer and in- 
ner edges) in Fig. 3(a) is 29. SA, while the width of the 
thinner ring in Fig 3(b) is 19. 9A. When observing the 
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FIG. 3: (color online) (a) Hexagonal zigzag quantum ring 
with Nout — 21 and Ni„ = 7. (b) Thinner hexagonal zigzag 
quantum ring, with the same Nout = 21, but with Nin ~ 12. 
(c) Energy spectrum as function of the magnetic flux for the 
structure in (a), (d) Energy spectrum as function of the 
magnetic flux for the structure in (b). 
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effects of varying the width in the energy spectra, these 
ring widths can be compared to the magnetic lengths for 
corresponding magnetic fluxes, as described in the pre- 
vious section. Corroborating the idea that the couphng 
between inner and outer edge states is directly related to 
the appearance of the energy gaps between sub-bands, 
we clearly see that the thinner quantum ring shows en- 
ergy gaps in the spectrum until higher values of magnetic 
fluxes (smaller magnetic lengths). 



B. Zigzag versus armchairm edges: differences in the 
quantum ring spectra around the Dirac point 

Having in mind the possible importance of the edge 
structure on the electronic structure of graphene quan- 
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tum rings, we now look to hexagonal quantum dot and 
quantum ring systems with inner and outer armchair 
edges (Figs. 4(a) and 4(b)). The corresponding elec- 
tronic structures as a function of magnetic fleld are shown 
in Figs. 4(c) and 4(d). The number of hexagonal pla- 
quetes in each side of the hexagonal dot considered is 
Nout = 13 (the counting for armchair edge terminations 
takes in account only the outermost plaquetes), corre- 
sponding to a total of 2814 atoms in the nanostructure. 
For the hexagonal ring, we considered N^ut = 13 and 
Nin = 6, where 7Vi„ is again the number of hexagonal 
cells in one side of the hexagon removed. 

Recalling that at i?=0 there are no states associated to 
armchair edges near the Dirac point [S], the central part 
(around E=Qi) of the quantum ring spectrum here is com- 
pletely different than in the case of zigzag edges (compare 
with Fig. 3). There are still edge states sub-bands de- 
fined, each one containing six energy levels, however the 
difference is that there is now a wide sub-band around 
the Dirac point. An interesting observation is the inter- 
change between electron-like and hole-like states in this 
region, as a function of magnetic field. The huge differ- 
ence in the electronic dispersion should be reflected in 
the related transport properties. 

Similarly to the zigzag case, a clear and strong local- 
ization of the charge density at the inner and outer edge 
occurs in the decoupled edges limit (high magnetic field), 
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FIG. 4: (a) Hexagonal graphene quantum dot with armchair 
edges, for which Nout ~ 13. (b) Hexagonal graphene quan- 
tum ring with armchair edges, with Nout = 13 and Nin = 6. 
(c) Energy spectrum for the structure in (a), (d) Energy 
spectrum for the structure in (b), where the arrows indicate 
states whose charge densities are plotted in Figure 5. 



FIG. 5: (color online) Electronic charge distribution of the 
selected edge states indicated by the arrows and correspond- 
ing letters in Fig. 4(d). (a) an outer edge state, (b) an inner 
edge state, (c) coupling between the inner and the outer 
edges, at an anticrossing. 
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as observed in the examples of Fig. 5(a) and Fig. 5(b), 
corresponding to the states pointed by the arrows (a) 
and (b) in the ring spectrum (Fig. 4(d)). Around an 
anticrossing, as for the state pointed by arrow (c), the 
charge density, observed in Fig. 5(c) is spread out on the 
two edges, indicating the edge coupling. As a difference 
between zigzag and armchair cases edges, we see that an 
armchair termination leads to a sublattice admixture of 
the charge density, different from the case of the zigzag 
edges, were there is a sublattice modulation [TU] . 



C. Asymmetries introduced by the corners in 
diamond rings 

Next we consider a rhombus-shaped (diamond) 
graphene quantum ring which has only zigzag edges. This 
ring is interesting because of its two-fold rotational sym- 
metry and because the upside outer edges of the dia- 
mond (in the perspective of the pictures in Figs. 6(b-d)) 
are from one sublattice meanwhile the downside outer 
edges are from the other sublattice. However, similarly 
to the triangular quantum ring, the up and down cor- 
ners belong to a different sublattice of the neighboring 
edges. The same sublattice effects occur at the inner 
edges. The number of atoms for this kind of ring is given 
by 2{Nl^^ - Nl) + 4(iV„„t - Ar,„) - 2, within our defini- 
tion. The energy spectrum of a diamond-like ring defined 
by Nout — 32 and iVi„ = 10 (1934 atoms) is plotted in 
Fig. 4(a). One can clearly observe the evolution of the 
two level bands, as expected from the symmetry of the 
structure. 

All the rings investigated in the present work with 
zigzag edges show similar energy spectra: at the low 
field limit sub-bands of energy states are formed, with 
the number of levels in each sub-band given by the ring 
symmetry and well defined outer and inner edge states at 
higher magnetic fields. Nevertheless, a closer look at the 
electronic charge density associated to different up and 
down going states, as well as at anticrossings, for these 
diamond-like rings reveals a further ingredient in the ef- 
fect of edges on the electronic properties of quantum dots 
and rings in graphene, namely the edge junctions at the 
corners. In a diamond-like ring the junctions between the 
zigzag edges define single armchair-like units at the left 
and right (inner and outer) edges (Figs. 6 (b)-(d)), while 
the upper and lower corners remain zigzag like. One can 
see, in the sequences of electronic charge distributions in 
Fig. 6, the high density around the armchair like corners, 
independently from being a state at the outer edge. Fig. 
6(b), inner edge. Fig. 6(c), or even at an anticrossing, 
Fig. 6(d). 

This situation calls the attention to the possible role 
of the edge junctions on the localization of the electronic 
charge in graphene nanostructures, i.e., the localization 
of the electronic charge at a rough interface may depend 
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FIG. 6: (color online) (a) Energy spectrum as function of the 
magnetic flux for a rhombus-shaped quantum ring with zigzag 
outer and inner edges, containing Nout = 32 and Nin = 10 
hexagonal plaquetes in each side length. Electronic charge 
densities: (b) going down state marked with the (b) arrow 
in the spectrum, (c) the going up state marked with the (c) 
arrow and (d) for an anticrossing state marked with the (d) 
arrow. The radii of the circles denote the magnitude of the 
charge density. 



on the symmetries at the corners that define the edge 
landscape. 



D. Round rings - effects of irregular edges 

We then analyze the cases of a round graphene dot 
and ring. Here the edges of the structures are irregu- 
lar and were defined in a way to best approach circular 
geometries for outer and inner edges, taking care not to 
leave edge atoms with only one nearest-neighbor [T3j, as 
observed in Figs. 7(a) and 7(b). Figures 7(c) and 7(d) 
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show the corresponding energy-magnetic flux spectra for 
these two structures. For the round dot, the number 
of atoms that has been taken is 2283, deflning a radius 
of « 47.1 A. For the round ring, the external radius is 
« 47.1 A, and the internal radius is « 7.3 A, containing 
a total of 2226 atoms. 

One can perceive from these spectra that, despite the 
irregularities of the edges, the main effects observed from 
the previous geometries are robust and keep present 
here. Comparing the spectra for the circular ring (Fig. 
7(d)) with the one for the circular dot (Fig. 7(c)), it 
is again clear that the circular antidot introduces inner 
edge states whose energies are increased with increasing 
magnetic flux. In the low flux limit of the ring spec- 
tra, anticrossing levels are again observed in the edge 
states region, indicating the coupling of inner and outer 
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FIG. 7: (a) Round graphene quantum dot with 2283 atoms, 
(b) Round graphene quantum ring with 2226 atoms, (c) En- 
ergy spectrum for the structure in (a), (d) Energy spectrum 
for the structure in (b). 



edge states, exactly as observed and described for the 
structures with well deflned edge structures (zigzag or 
armchair). The main difference is that, as for this geom- 
etry there is no rotational symmetry, we do not observe 
the formation of edge-states sub-bands with well defined 
number of energy levels. 

V. CONCLUSIONS 

The present paper focus on the single particle elec- 
tronic properties of finite graphene structures. The be- 
havior of edge states in graphene rings is investigated, 
through the numerical calculation of the electronic en- 
ergy spectra of these rings as a function of a perpen- 
dicular magnetic field and the mapping of charge den- 
sity distributions. Several similar patterns may be found 
among quantum rings with different symmetries (trian- 
gular, hexagonal, diamond shapes), including the forma- 
tion of sub-bands of edge states energy levels, separated 
by energy gaps (anticrossings) . The choice of quantum 
rings revealed a strategic one because of the clear relation 
between the symmetry of the structure and the number 
of levels in the edge states sub-bands. Furthermore, the 
edge states levels within the sub-bands can be perfectly 
associated to either inner or outer edges, as well as the 
"bulk" region of the structure (couphng between edges) , 
therefore enabling a good framework for studying the in- 
fluence of the edges on the electronic structure and charge 
distribution. If edge terminations (zigzag or armchair) 
show to play an important role on the electronic prop- 
erties, specially for the states around the Dirac point, 
the junction of the edges (corners) can also be crucial for 
charge density localization patterns. 
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